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Abstract 

The method of the large mass expansion (LME) has the technical advantage 
that two-loop integrals occur only as bubbles with large masses. In many cases only 
one large mass occurs. In such cases these integrals are expressible in terms of V- 
functions, i.e. they can be handled completely analytically avoiding even recursions 
and therefore this approach may find a wide field of application. We consider it 
necessary to investigate the precision of this method and test it for several two- 
loop vertex functions ocurring in the Z — > bb decay by comparing it with the small 
momentum expansion. It turns out that in general high order approximants have 
to be taken into account for a sufficient accuracy. 
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1 Introduction 



The experiments at LEP/CERN and SLC/SLAC have reached a precision, which goes 
beyond all former expectations. It turns out that we are approaching the limits of our 
theoretical understanding and in order to fully evaluate the present and future experimen- 
tal data (LEP1 data still to be investigated, LHC and a possible e + e~ linear collider) a 
detailed analysis of higher order corrections is necessary. In particular at higher energies 
radiative corrections become more and more important since they can grow strongly. 

Quite a number of groups have started to develop methods for the evaluation of two- 
loop diagrams. Analytic results can in general only be obtained for specific mass relations, 
like e.g. all masses zero or only one non-zero mass in a diagram. Since in the two-loop 
order for most of the interesting processes the number of diagrams can become of the order 
of 1000, it is also necessary to calculate these diagrams with extremely high precision for 
the many different masses occurring in the Standard Model (SM). 

A method which provides particular high precision for many kinematical configurations 
of interest is that of Taylor expansion in an external momentum squared improved by 
conformal mapping and the Pade summation technique Another method, which is 
particularly convenient to program in a formulae manipulating language (we use FORM 
0), is the large mass expansion (LME) |5). If we have, e.g., diagrams with one large 
mass like a top-quark in a virtual line j|, [5], the method is shown to be appropriate. To 
obtain a satisfying precision, however, it turns out to be absolutely necessary that higher 
order terms of the LME are taken into account. For two- loop propagator- type diagrams 
such an analysis was carried out in ||. In this paper we want to consider three-point 
functions. The particular process we have in mind is Z — > bb, i.e. we investigate vertex 
diagrams contributing to this process which contain a top-quark and compare the results 
of the above two methods in order to find out to what precision it is possible to obtain 
results from the LME. In most of the cases under consideration the Taylor expansion 
needs only up to 9 coefficients for a very high precision since we can confine ourselves 
here to momenta squared below the lowes threshold. 

2 Large Mass Expansion (LME) 

Two-loop diagrams as shown in Fig.l with two different masses on virtual lines, one of 
which a top, arise in the process Z — > bb. W and Z are the gauge bosons with masses M w 
and M z , respectively; <fi is the charged would-be Goldstone boson (we use the Feynman 
gauge); t and b are the t- and b-quarks. In our analysis we neglect the 6-quark mass 
(rrib = 0). Thus we have the following kinematics (see Fig.l): p\ = p\ = 0, (p2 — Pi) 2 = <? 2 
and q 2 = M§ for the on-shell Z. For such type of diagrams analytic results are impossible 
to obtain with presentday technology. 

Applying the method of the LME, My/ and Mz are considered as small. In the 
framework of this expansion contributions from additional subgraphs are to be taken 
into account together with the Taylor expansion of the initial diagrams w.r.t. external 
momenta and light masses. These subgraphs restore the analytic properties of the initial 
diagrams (like logarithmic behaviour). In other words, the Taylor expansion of the initial 
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diagrams produces extra infrared singularities which are compensated by singularities of 
the additional subgraphs so that only the singularities of the original diagrams survive. 
For a given scalar graph T the expansion in the large mass is given by the formula 

F r (q, M, m, e) M ~°° ^ F V / 7 (q, m, e) o T q ^ m -,F 1 {q 1 , M, m 7 , e), (1) 

7 

where 7's are subgraphs involved in the LME, r/7 denotes shrinking of 7 to a point; F 1 is 
the Feynman integrand corresponding to 7; T qj)mj is the Taylor operator expanding this 
integrand in small masses {m 7 } and external momenta {g 7 } of the subgraph 7 ; o stands 
for the convolution of the subgraph expansion with the integrand -Pr/ 7 - The sum goes 
over all subgraphs 7 which (a) contain all lines with large masses, and (b) are one-particle 
irreducible w.r.t. light lines. 

At this point it becomes clear what the difference is between the small-g 2 expansion 
and what is called here the LME: in the former case we assume all masses large, i.e. 
q 2 <C M^Mfj-mj 3 while in the latter case only m t is considered as large and all other 
parameters small, i.e. q 2 , M^, M\ <C m 2 . In this sense both methods are LM expansions. 
The technical advantage of the second method is, however, that only bubbles with one 
mass occur, which can be expressed in terms of T functions, while in the case of the small 
q 2 expansion bubbles with different masses are involved, which are much more difficult 
to evaluate. Of course also the number and structure of the subgraphs is different in the 
two cases. 

The structure of the LME of the diagrams under consideration is given in Fig.2. Bold, 
thin and dashed lines correspond to heavy-mass, light-mass, and massless propagators, 
respectively. Dotted lines indicate the lines omitted in the original graph T to yield the 
subgraph 7. r/7 (see (|I|)) then consists out of all the dotted lines after schrinking 7 
to a point. Thus, in case the last and in case N the two last contributions vanish in 
dimensional regularization. Finally the LME of the above diagrams has the following 
general form: 

1 N n / M 2 y / 2 \ 3 m 2 

^ = AE £ R Pi E Ai,*(9 2 , M 2 W , » 2 ) ln k ^ (2) 

where m is the highest degree of divergence (ultraviolet, infrared, collinear) in the various 
contributions to the LME (m < 3 in the cases considered). M^/m 2 and q 2 /m 2 are consid- 
ered as small parameters. Aij t k are i n general complicated functions of the arguments, i.e. 
they may contain logarithms and higher polylogarithms (see also the explicit examples 
below). (0) implies that the difference F — F^ of the integral F and its approximation 
behaves like 

F-FZ = 0[ K W \ K \i +2 * , withi+j = N. (3) 

V ( m t) J 

One might expect this difference to indicate the order of the real error of the asymptotic 
expansion. For q 2 small this will indeed be the case, for q 2 ~ Mf the "small" parameters 
are not really small anymore and the validity of the LME can only be inferred from 
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the comparison with results from other methods, for which we take here the small q 2 
expansion. In particular nothing can be said about a range of convergence of the LME. 
We can, however, at least formally apply the Pade approximation method to improve the 
convergence, the encouraging results of which will be demonstrated below. 

3 Results 

The small momentum expansion of cases 7, 7.1 and 7.2 (notation due to FST of Q) is 
described in detail in Ref . [0] . The additional subgraph arising in these cases is the same 
irrespectively of the mass distribution and is shown in Fig. 2 (case 7). 

In the case of the LME two additional subgraphs arise in each of the cases 7.1 and 
7.2. Beyond that these sets of additional subgraphs are also of quite different structure. 
Furthermore the 2 nd and 4 th graphs in the row (see Fig. 2) produce 1/e 3 terms which 
cancel, however. Since there are no UV divergences, these must be mixtures of infrared 
and collinear ones. They determine the highest degree of divergence in these cases and 
thus the highest power of the logarithm as discussed in (2). 

Our numerical results for cases 7.1 and 7.2 are presented in Fig.3, and for q 2 = M| in 
Table 1. In the figures we show the small g 2 -expansion in comparison with the lowest order 
approximation and the sum of terms with N=9, see (2). The scale parameter fi = m t , i.e. 
only k=0 contributes in (2). We see that up to q 2 = Mf the sum of 11 terms agrees quite 
well with the small g 2 -expansion, while for higher q 2 the agreement very quickly worsens. 
For this reason we formally apply also the Pade summation technique. With 11 terms in 
the series, a [5/5]-approximant (long-dashed) can be constructed. It is seen that indeed 
this improves the situation considerably up to the first threshold though for a better 
agreement many more terms in the LME would be needed. In the small g 2 -expansion 
only 9 terms were taken into account, i.e. a [4/4] approximant is calculated. The high 
precision of the results in this case even up to the threshold (in case 7.2) is due to the 
mapping applied in addition (see [|]]). From Table 1 we see that for q 2 = Mf indeed a 
rather precise result can be obtained with 11 terms from the LME, in particular if Pade's 
are applied. 

In case an extremely high agreement between the small-g 2 expansion and the LME 
is obtained (see Table 1). No deviation in a figure like above would be seen. 

The last two cases: A (for asymmetric) and N (for nonplanar) are difficult from the 
point of view of the small g 2 -expansion. While for the planar and symmetric diagrams 
a "standard" numerical program was written (applied e.g. in 0), this program is not 
without considerable changes applicable for an asymmetric diagram and not without even 
more extensions for a nonplanar one. Thus in these cases just the analytic expressions 
for the first 9 Taylor coefficients of the small g 2 -expansion were produced with FORM 
in terms of bubble integrals and then numerically evaluated. The problem is that higher 
coefficients get quite lengthy in this manner so that the FORTRAN programs are difficult 
even to compile (in the multiple precision version of D. H. Bailey ||)! 

In case N the lowest threshold is indeed quite high and therefore the obtained coeffi- 
cients yield an extremely precise result for the diagram at q 2 = Mf . What concerns the 
LME, the situation is essentially the same as in cases 7.1 and 7.2 (see Fig.3 and Table 1). 



3 



Case A is interesting due to the fact that the lowest threshold is rig ht at q 2 = Mf m 
fact the singular contribution at this point should cancel against some Bremsstrahlung 
contribution to the Z — > bb decay). Table 2 demonstrates that in this case the LME 
yields results of extremely high precision, which perfectly agree with the results from the 
small-g 2 expansion at low q 2 . Near and above the threshold, however, the precision of 
the small g 2 -expansion with a total of 9 coefficients is getting worse, but nevertheless the 
agreement is still quite surprising. In this case, obviously, the LME is superior to the 
small momentum expansion. 

_p(o) j n Tables 1 and 2 corresponds to the lowest order of the diagrams (N=-l or 0, 
respectively) and to the order of expansion performed in Ref. ||. 



Table 1: Values of diagrams at q 2 = 90 2 , ji = m t = 180, Mw = 80. 
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Table 2: Case A for different q 2 . 
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A Analytical Expressions 

In this section we present the first two terms of the LME for the diagrams under consid- 
eration. All Feynman integrals are normalised as follows 

(4tt) 2£ T 2 (l + e) T 2 (l - e) 4e f d d h d d k 2 1 



(16tt 2 ) 2 M 4 r(l-2e) ^ J i{2n) d i{2Ti) d DxD 2 D 3 D A D 5 D 

where Di = (q 2 — m 2 ) are corresponding inverse propagators and d = 4 — 2e is the space- 
time dimension. Instead of the large mass M, the small mass m and the momentum 
squared, q 2 , we introduce the dimensionless parameters z = m 2 /M 2 and s = q 2 /m 2 and 
write the results of the LME as a series in z with each coefficient being s-dependent. 

^7.1 = ~ ~ -{-Lg ~ l m + 1) - ~{L g + L M f + (L q + L M ) - 1 



(^2 - ~ L i + 2 L ?) ( Lm ~ Lm + | o 
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+ (i-L,)(i(Li,-L 



'2^m + 



24 + s 
24 



2 9 



24 + s 8 - 3s 

——L M + L m + ^- 

7s - 16 16 - 23s 

f" -^m H o H 77 

8 16 



+ 0(z 2 ), (5) 



zs J 7 . 2 = 



{^-- £ L q + lLl)(L M -L m + h(s)-i) 

(I - L q ) (h 2 m - \h M + L m - h{s)L m + J 2 (s) -l)-h L ' M - L* m ) 



-2 m 2 

r t^m - - ^i(s)^ + / 2 (s)L m - I 3 (s) + 1 + O(z) 



1 

'2" 



J = -l + C(2)-^ 



L m + 



144 + 49s 4 + s 
72 2~ 



C(2) 



+ 0(^ 2 ), 



(6) 
(7) 



zsJ A = +y £ log ^ ~ s ) " 2 (^ Lm + log ^ ~ S ^ + l0g2 ^ ~ ^ + Li2 ^ S 



+ 



24e 
1 + s 
24 



L((i + a )io g (i_ a ) 



$ $ $ 

-L m + ^L M + - 



1 3 

J7V ""27"2 +jLm_2; 



L m + L M - 1) log(l - s) + log 2 (l - s) + Li 2 (s 
- {Lm — Lm + 



72 + 7s x 1 2 

~2 m _ l^rn^M 



3 2 



144 + 7s 
24 



Lm + 



48 

720 + 181s 



288 



+ 2C(2) 



+ 0(z 2 ), 



(9) 



where Lm = log(M 2 //i 2 ), L m = log(m 2 / ' ji 2 ) and L q = log(— g 2 /^ 2 ) with /i being an 
arbitrary scale parameter. The functions I n (s) can be expressed in terms of polylogarithms 
Lifc. Here we give only their integral representations 



In(s) 



1 ^ 

2~d~ 



log"(l -ts/4) 



dt. 



(10) 
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Figure 1: Two-loop diagrams with two different masses in internal lines arising in the 
process Z — ► 66. 
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Figure 2: The structure of the LME, see explanations in the text. 
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Figure 3: Results of the LME for the real finite parts of diagrams 7.1, 7.2 and N with 
/i 2 = m 2 . Solid curves represent the small-g 2 expansions, middle-dashed the leading term 
of the LME, short-dashed the sum of 11 terms in the LME, long-dashed the [5/5] Pade 
approximant from the LME. 




